On the problem of inflation in nonlinear multidimensional cosmological models 
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We consider a multidimensional cosmological model with nonlinear quadratic and quartic 
actions. As a matter source, we include a monopole form field, D-dimensional bare cosmological 
constant and tensions of branes located in fixed points. In the spirit of the Universal Extra Dimen- 
sions models, the Standard Model fields are not localized on branes but can move in the bulk. We 
define conditions which ensure the stable compactification of the internal space in zero minimum of 
the effective potentials. Such effective potentials may have rather complicated form with a number 
of local minima, maxima and saddle points. Then, we investigate inflation in these models. It is 
shown that E? and 7?* models can have up to 10 and 22 e-foldings, respectively. These values are not 
sufficient to solve the homogeneity and isotropy problem but big enough to explain the recent CMB 
data. Additionally, i?* model can provide conditions for eternal topological inflation. However, the 
main drawback of the given inflationary models consists in a value of spectral index Us which is less 
than observable now Us ~ 1. For example, in the case of 7i* model we find Ug ~ 0.61. 

PACS numbers: 04.50.+h, 11.25.Mj, 98.80.-k 



INTRODUCTION 



Recently, very elegant idea of inflation has achieved 
f!^ , spectacular success in explaining the acoustic peak 

Q ■ structure seen in CMB (see e.g. Jj]). It is very dif- 

(— I I ficult to explain correctly the observable large-scale 

structure formation without taking into account the 
stage of early inflation. There is a big number of dif- 
' ferent models of inflation. However, the most of them 

^ ' are poorly related to fundamental physics. In these 

^ : models, the stage of inflation occurs due to a spe- 

cial form of scalar field potential. Here, the origin of 
' scalar field and the form of its potential is usually re- 

, mained out of the scope of these papers. However, 

^\ • it is well known that scalar fields has naturel origin 

' from higher-dimensional theories. They are geomet- 

OO ! rical moduli (radions, gravexcitons) which are related 

. to the shape of internal spaces (to scale factors of the 

internal spaces). After dimensional reduction to four 
dimensions, scalar field potential is completely defined 
, by the topology and matter content of original higher- 

<^ • dimensional model 0, 0]. Thus, it is of undoubted 

' interest to realize inflation in these models (see e.g. 

[1, 0, [6^ in string theory and in multidimensional 
cosmological models and references therein). 

On the other hand, scalar fields with corresponding 
potentials originate naturally from nonlinear gravita- 
tional models where Lagrangian is a function of scalar 
curvature f{R)- It is well known that such models are 
equivalent to linear ones plus scalar field. This scalar 
field corresponds to an additional degree of freedom 
of nonlinear models. For motivation of these theories, 
see review [8]. For example, among others higher- 
order gravity theories, f{R) theories are free of ghosts 
and of Ostrogradski instability [^. These models at- 
tract great attention because can provide the late-time 



acceleration of our Universe due to a special form of 
scalar field potentials (see e.g. d, [13, El and refer- 
ences therein), which is interesting alternative to the 
cosmological constant. These models can also provide 
the stage of early inflation both in four-dimensional 
(see the pioneering paper by Starobinsky jl2] and nu- 
merous references in u^) and multidimensional 
[ni,[il[il cases. 

In our paper we combine both of these approaches. 
Namely, we consider multidimensional models with 
nonlinear action. To start with, we investigate the 
most simple linear multidimensional model. We show 
that such model can provide power-law inflation. Un- 
fortunately, it takes place for the branch where in- 
ternal space is decompactified ^. Then, to get infla- 
tion of the external space with subsequent stabiliza- 
tion of the internal spaces, we turn to multidimen- 
sional nonlinear models with quadratic and quartic 
nonlinearities. First, we obtain the conditions of the 
internal space compactification (stabilization). Sec- 
ond, for corresponding effective potentials, we investi- 
gate the possibility of the external space inflation. We 
show that in the quadratic and quartic models we can 
achieve 10 and 22 e-folds, respectively. These num- 
bers are sufficient to explain the present day CMB 
date, but not enough to solve the horizon and flat- 
ness problems. However, 22 e-foldings is rather big 
number to encourage the following investigation of 
the nonlinear multidimensional models to find theo- 
ries where this number will approach to 50-60 e-folds. 
Even more, this number (50-60) can be reduced in 
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Each time when we consider multidimensional cosmological 
models we should remember about the problem of the inter- 
nal space stabilization/compactification. If such stabilization 
is absent we confront with the variation of four-dimensional 
fundamental constants. General method of the internal space 
stabilization was described in Q and was applied after that 
to numerous models. In the present paper, we follow also this 
method. 
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models with long matter dominated stage followed by 
inflation with subsequent decay into radiation. Pre- 
cisely this scenario takes place for our models where 
we find that e-folds can be reduced by 6 if the mass 
of decaying scalar field m ~ ITeV. So, we believe that 
the number of e-folds is not a big problem for pro- 
posed models. The main problem consists in spectral 
index ~ 0.61 (for the quartic model) which is less 
than observable « 1. A possible solution of this 
problem may consist in more general form of the non- 
linearity f{R)- For example, it was observed in [T5| 
that simultaneous consideration quadratic and quar- 
tic nonlinearities can flatten the effective potential. 
We postpone the investigation of this problem for our 
following paper. 

To conclude, we want to indicate two interesting 
features of models under consideration. First, the 
quartic model can provide the topological inflation. 
Here, due to quantum fluctuation of scalar fields, in- 
flating domain wall has fractal structure (inflating do- 
main wall will contain a number of new inflating do- 
main walls and each such domain walls will contain 
again a new inflating walls etc [16]). So, we arrive at 
the eternal inflation. Second, obtained solution has 
the property of the self-similarity transformation (see 
Appendix B). It means that in the case of zero mini- 
mum of the effective potential and fixed positions for 
extrema in (t/j, 0)— plane, the change of the hight of ex- 
trema results in rescaling of the dynamical character- 
istics of the model (graphics of the number of e-folds, 
scalar fields, the Hubble parameter and the acceler- 
ation parameter versus synchronous time) along the 
time axis. The decrease (increase) of hight in c times 
(c is a constant) leads to the stretch (shrink) of these 
figures along the time axis in -^/c times. 

The paper is structured as follows. In Sec. II we 
consider the linear (on scalar curvature) model. The 
nonlinear quadratic E? and quartic i?^ models are in- 
vestigated in Sec. Ill and Sec. lY, respectively. Here, 
we find the range of parameters where the internal 
space is stabilized and investigate a possibility for the 
external space inflation. A brief discussion of the ob- 
tained results is presented in the concluding Sec. Y. 
The Friedmann equations for multi-component scalar 
field models are reduced to the system of dimension- 
less first order ODEs in Appendix A. In Appendix 
B, we show that the dynamical characteristics (e.g. 
the Hubble parameter, the acceleration parameter) of 
considered nonlinear models satisfy the self-similarity 
condition. 



II. LINEAR MODEL 



space-time (usually, we have in mind that Dq ~ 4) 
and Ml corresponds to di -dimensional internal space 
which is a flat orbifold^ with branes in fixed points. 
Scale factor of the internal space depends on coordi- 
nates X of the external space-time: ai(x) — Lpie^ 
where Lpi is the Planck length. 

First, we consider the linear model f{R) = R with 
/^-dimensional action of the form 



S 



1 



d''x^\^\[R[g^''^]-2KD]+Sm+Sb. 



(2.2) 

Kd is a bare cosmological constant'^. In the spirit 
of Universal Extra Dimension models [l3], the Stan- 
dard Model fields are not localized on the branes but 
can move in the bulk. The compactification of the 
extra dimensions on orbifolds has a number of very 
interesting and useful properties, e.g. breaking (su- 
per)symmetry and obtaining chiral fcrmions in four 
dimensions (see e.g. paper by H.-C. Cheng at al in 
[T3|). The latter property gives a possibility to avoid 
famous no-go theorem of KK models (see e.g. [l8j|). 
Additional arguments in favor of UED models are 
listed in [19]. 

Following a generalized Freund-Rubin ansatz [23] 
to achieve a spontaneous compactification M M = 
Mq X Ml, we endow the extra dimensions with real- 
valued solitonic form field i^^^^ with an action: 

s.^4//°xv^5^(f">)'- (") 

This form field is nested in di-dimensional factor space 
Ml, i.e. -F'^' is proportional to the world- volume 
of the internal space. In this case (l/di!) (^^^^) = 

fi ja^i^ , where /i is a constant of integration [13] . 

Branes in fixed points contribute in action func- 
tional (|2?2|1 in the form Qf^: 



Sh — 



E 

fixed 



Mo 



(2.4) 



where g^^'^lx) is induced metric (which for our geom- 
etry (12. 1|) coincides with the metric of the external 
space-time in the Brans-Dicke frame) and Lf, is the 
matter Lagrangian on the brane. In what follows, we 
consider the case where branes are only characterized 
by their tensions ib(fc) — —T(k) : k = 1,2, ... ,m and 
m is the number of branes. 



To start with, let us define the topology of our mod- 
els. We consider a factorizable ZJ-dimensional metric 



g(^)=g(°)(.)+il,,e^^^(-)gW, (2.1) 

which is defined on a warped product manifold M — 
Mq X Ml. Mq describes external Do-dimensional 



^ For example, jZi and T'^ /Z2 which represent circle and 
square folded onto themselves due to Z2 symmetry. 

^ Such cosmological constant can originate from D-dimensional 
form field which is proportional to the D-dimensional world- 
volume: F^^^-O = (C/y|g(-D)|)£MiV...Q^ this case the 
equations of motion gives C = const and term in action 
is reduced to {1/ D\)Fmn ...qF^'"^ ■'■'^ = -C^- 
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Let Pq be the internal space scale factor at the 
present time and f3^ = (3^ — /3q describes fluctuations 
around this value. Then, after dimensional reduction 
of the action (|2.ip and conformal transformation to 



-2/(Do-2) 



the Einstein frame g'jj^J = (^e'''^^ ^ 
arrive at effective Z?o-dimensional action of the form 



~(o) 

glJ, we 



1 

2«^0 J Mo 



g^'^^'^^d^ipd.ip~2U,ffiip)} , (2.5) 



where scalar field ip is defined by the fluctuations of 
the internal space scale factor: 



^ V ^0-2 ^ 



(2.6) 



and G :— Kq/Stt := n'jj / {8nVdi) {Vd^ is the internal 
space volume at the present time) denotes the Dp- 
dimensional gravitational constant. The effective po- 
tential Ueff{(p) reads (hereafter we put Dq — 4): 



/ 2di 
V dl+2 



A 



D 



fie 



V 



Ae 



(2.7) 



where = K|)/iVa(o)\ and A = -k^ X^feLi ^k)- 

Now, we should investigate this potential from the 
point of the external space inflation and the inter- 
nal space stabilization. First, we consider the latter 
problem. It is clear that internal space is stabilized if 
Ueffi^) has a minimum with respect to f. The po- 
sition of minimum should correspond to the present 
day value if = 0. Additionally, we can demand that 
the value of the effective potential in the minimum 
position is equal to the present day dark energy value 
Ueff{(p = 0) ~ Ade ^ 10~^''cm~^. However, it re- 
sults in very flat minimum of the effective potential 
which in fact destabilizes the internal space [2^ . To 
avoid this problem, we shall consider the case of zero 
minimum Ueffi^p = 0) = 0. 

The extremum condition dUeff /dip\^^g = and 
zero minimum condition Ueff^^p = 0) = result in 
a system of equations for parameters Kd, fi and A 
which has the following solution: 



Ad = fl = A/2 . 



(2.8) 



For the mass of scalar field excitations (gravexci- 
tons/radions) we obtain: w? = d^Uef f / dip'^\^_^ ~ 

{Mi/{di + 2)) Kd In Fig. 1 we present the effec- 
tive potential (|2.7p in the case di = 3 and Au = 10. 
It is worth of noting that usually scalar fields in the 



present paper are dimensionless^ and f/e//, A^, , A 
are measured in Mp^ units. 

Let us turn now to the problem of the external space 
inflation. As far as the external space corresponds to 
our Universe, we take metric 5^"^ in the spatially flat 
Friedmann-Robertson- Walker form with scale factor 
a{t). Scalar field Lp depends also only on the syn- 
chronous/cosmic time t (in the Einstein frame). 

It can be easily seen that for 93 >> (more pre- 
cisely, for ip > (firnax = \/ {di + 2) /2di In 3) the poten- 
tial (|2.7p behaves as 



[/e//((p)«Az3e-v^^, 



with 



q ■= 



2di 
di+2 



(2.9) 



(2.10) 



It is well known (see e.g. [H, [H, HI, HI]) that for 
such exponential potential scale factor has the follow- 
ing asymptotic form: 



a{t) - t^/g 



(2.11) 



Thus, the Universe undergoes the power-law inflation 
a q < 2. Precisely this condition holds for eq. (|2.10p 
if di > 1. 

It can be easily verified that (p > (pmax is the 
only region of the effective potential where inflation 
takes place. Indeed, in the region ip < the lead- 
ing exponents are too large, i.e. the potential is too 
steep. The local maximum of the effective potential 
(4/27)Ai3 at 



U, 



= y/{di+2)/2di In 3 
is also too steep for inflation because the slow-roll pa- 



rameter 77,, 



_ 3rfi 



ii+2 



^ 1 < 



rjmax I < 3 and does not satisfy the inflation condition 
77I < 1. Topological inflation is also absent here be- 
cause the distance between global minimum and local 
maximum Lpmax — V(di + 2)/2di In 3 < 1.35 is less 
than critical value ipcr > 1-65 (see [H, [13, E^)- It 
is worth of noting that rjmax and ipmax depend only 
on the number of dimensions di of the internal space 
and do not depend on the hight of the local maximum 
(which is proportional to A^i). 

Therefore, we have two distinctive regions in this 
model. In the first region, at the left of the maxi- 
mum in the vicinity of the minimum, scalar field un- 
dergoes the damped oscillations. These oscillations 
have the form of massive scalar fields in our Universe 
(in [3] these excitations were called gravitational ex- 
citons and later (see e.g. (Hj) these geometrical mod- 
uli oscillations were also named radions). Their life- 
time with respect to the decay (/s — > 27 into radia- 
tion is [13, HE 111 r ~ {Mpi/mfTpi. For example. 



* To restore dimension of scalar fields we should multiply their 
dimensionless values by Mpi/^/Stt. 
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FIG. 1: The form of the effective potential (|2^ in the case di = 3 and Ad = /f = A/2 = 10. 



we obtain r ~ 10 s, lO'^s for m ~ lOTeV, lO^TcV 
correspondingly. We remind that in our case = 
{Adi/{di + 2))Ad- Therefore, this is the graceful exit 
region. Here, the internal space scale factor, after de- 
cay its oscillations into radiation, is stabilized at the 
present day value and the effective potential vanishes 
due to zero minimum. In second region, at the right 
of the maximum of the potential, our Universe under- 
goes the power-low inflation. However, it is impossible 
to transit from the region of inflation to the graceful 
exit region because given inflationary solution satisfles 
the following condition (p > 0. There is also serious 
additional problem connected with obtained inflation- 
ary solution. The point is that for the exponential 
potential of the form (|2.9p . the spectral index reads 
[23,i25,]5: 



2-3q 

2-g 



(2.12) 



In our case (|2.10p . it results in ris = 1 — di- Obviously, 
for di > 1 this value is very far from observable data 
Us 1. Therefore, it is necessary to generalize our 
linear model. 



2.2|1 ^. Then, effective potential takes the form 

Ueff {if A) = e"^^'^ 

/ 2di 



(3.1) 



where we put Ajj = in 

However, it is well known that scalar field (j) can 
naturally originate from the nonlinearity of higher- 
dimensional models where the Hilbert-Einstein linear 
lagrangian R is replaced by nonlinear one f{R). These 
nonlinear theories are equivalent to the linear ones 
with a minimal scalar field (which represents addi- 
tional degree of freedom of the original nonlinear the- 
ory). It is not difficult to verify (see e.g. [il,[2l|) that 
nonlinear model 



1 



1 

2 ./M 



(3.2) 



III. NONLINEAR QUADRATIC MODEL 

As follows from the previous section, we want to 
generalize the effective potential making it more com- 
plicated and having more reach structure. Introduc- 
tion of an additional minimal scalar field is one of 
possible ways. We can do it "by hand" , inserting mini- 
mal scalar field (j) with a potential U {(p) in linear action 



is equivalent to a linear one with conformally related 
metric 



(D) 
9ab 



2A4,/{D-2)-{D) 



(3.3) 



plus minimal scalar field cj) — hi[df /dR]/A with a po- 
tential 



1 



U{^) - ^e-^"^ [R{^)e^'^ - f{Rm] , 



(3.4) 



^ With respect to conformal time, solution 1 12.111 reads a{r]) ~ 
where /3 = —{4 — q)/(2 — q). It was shown in [33l | that for 
such inflationary solution (with q < 2) the spectral index of 
density perturbation is given by = 2/3 -|- 5 resulting again 
in (|2l^ . 



^ If such scalar field is the only matter field in these models, 
it is known (see e.g. fA Il3|) that the effective potential can 
has only negative minimum, i.e. the models are asymptoti- 
cal AdS. To uplift this minimum to nonnegative values, it is 
necessary to add form-fields [2l| . 
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where A = {D ~ 2) / {D - 1) = ^ {di + 2) / di + 'i 

and B = D/^{D^2){D-l) = A{di + 4)/{di + 2). 
After dimensional reduction of this Hnear model, we 
obtain an effective Dg-dimensional action of the form 



1 



(3.5) 



with effective potential exactly of the form p.ip . It is 
worth to note that we suppose that matter fields are 
coupled to the metric g'^^-' of the linear theory (see also 
analogous approach in [33|)- Because in all considered 
below models both fields ip and are stabilized in the 
minimum of the effective potential, such convention 
results in a simple redefinition/rescaling of the mat- 
ter fields and effective four-dimensional fundamental 
constants. After such stabilization, the Einstein and 
Brans-Dicke frames are equivalent each other (met- 
rics g^'^^ and g^'^^ coincide with each other), and linear 
g^^^ and nonlinear g*-^-* metrics in (I3.3P are related 
via constant prefactor (models became asymptotically 
linear)^. 

Let us consider first the quadratic theory 



f{R) = R + iR^ ~ 2Ai 



(3.6) 



For this model the scalar field potential (|3.4p reads 



2A7 



(3.7) 



It was proven [7] that the internal space is stabi- 
lized if the effective potential ()3.1|) has a minimum 
with respect to both fields (p and (j). It can be easily 
seen from the form of Uef fif-, (j}) that minimum (/)o 
of the potential U{(j)) coincides with the minimum of 
U,ff{ip,cj,) : dU/d^\^^ = ^ 90C/e//|^„ = 0. For 
minimum U{(j)Q) we obtain [l3l |: 



uiM = [{xo - 1) 

where we denote the constant xq 

(a 



-mi 



(3.8) 



= exp(yl0o) 
/{2A-B). 



^ B + y/W+(2A^^B)B8^D j / {2A- B). It 

is the global minimum and the only extremum of 
U{(j)). Nonnegative minimum of the effective po- 
tential Ueff takes place for positive £,,Ad > 0. If 
^,Ad > 0, the potential U{(j)) has asymptotic behav- 
ior U{(j)) +00 for (j) ±00. 

The relations (|2.8p , where we should make the sub- 
stitution Ad ^ (j (0o), are the necessary and sufficient 
conditions of the zero minimum of the effective poten- 
tial Ueff {if, (f)) at the point (ip = 0,(f) = ^o). Thus, 



if parameters of the quadratic models satisfy the con- 
ditions U{(j)o) = ff = A/2, we arrive at zero global 
minimum: Ueff{0,4>o) = 0. 

It is clear that profile = of the effective poten- 
tial Ueff has a local maximum in the region > 
because Ueffi^p, = <po) ^ if 1^9 ^ -f-oo. Such pro- 
file has the form shown in Fig.l. Thus, the effective 
potential Ueff has a saddle point {(p — Pmax, 4> = 0o) 
where ipmax — \/ {d\ + 2)/2di In 3. At this point 
Uefs\^ax = (4/27)[/(0o). The Figure [21 demonstrates 
the typical contour plot of the effecti ve p otential (j3.ip 
with the potential [/((/>) of the form (|3.7|) in the vicin- 
ity of the global minimum and the saddle point. 

Let us discuss now a possibility of the external space 
inflation in this model. It can be easily realized that 
for all models of the form (|3.1|) in the case of local 
zero minimum at (1^ = 0, 0o)i the effective potential 
will also have a saddle point at (1^9 = 'p>max, (t>o) with 
V(di +2)/2di In 3 < ^cr = 1-65 and the 
slow-roll parameter \rj^\ in this point cannot be less 
than 1: \'q^\ = 3di/{di -|- 2) > 1. Therefore, such sad- 
dles are too steep (in the section <j> = (pg) for the slow- 
roll and topological inflations However, as we shall see 
below, a short period of De Sitter-like inflation is pos- 
sible if we start not precisely at the saddle point but 
first move in the vicinity of the saddle along the line 
P ~ 'pmax with subsequent turn into zero minimum 
along the line (j) (pQ. Similar situation happens for 
trajectories from different regions of the effective po- 
tential which can reach this saddle and spend here a 
some time (moving along the line p « Pmax)- 

Let us consider now regions where the following con- 
ditions take place: 



U{^) » ffe 



(3.9) 



For the potential (|3.7p these regions exist both for 
negative and positive (p. In the case of positive (j) with 
exp(A0) > max {l , (B^Ai))^/^} we obtain 



(3.10) 



However, small quantum fluctuations around the minimum 
of the effective potential distinguish these metrics. 



where q is defined by Eq. ([2T0l) . qi := {2 A ~ Bf = 
d\/[{di + 2){di + 3)] and q> qi. For potential (fHTTOl) 
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FIG. 2: Contour plot of the effective potential U^f f{ip,(i>) p.ip with potential [/(<;/>) of the form (|3.7p for parameters 
di = 1 ,i^A_D = 1 and relations U{4>o) ~ fi = A/2 . This plot clearly shows the global minimum and the saddle. The 
colored lines describe trajectories for scalar fields starting at different initial conditions. 



the slow-roll parameters are^: 

e« ?7i « 772 « - + y (3.11) 

and satisfy the slow-roll conditions 6,771,7/2 < 1. As 
far as we know, there are no analytic solutions for such 
two-scalar-field potential. Anyway, from the form of 
the potential (|3.10[) and condition q > qi we can get an 
estimate a ^ t" with s >2/q (e.g. 2/q = 3,2,5/3 for 
di = 1,2,3, respectively). Thus, in these regions we 
can get a period of power-law inflation. In spite of a 
rude character of these estimates, we shall see below 
that external space scale factors undergo power-law 
inflation for trajectories passing through these regions. 

Now, we investigate dynamical behavior of scalar 
fields and the external space scale factor in more de- 
tail. There are no analytic solutions for considered 
model. So, we use numerical calculations. To do it, 
we apply a Mathematica package proposed in ad- 
justing it to our models and notations (see Appendix 
A). 

The colored lines on the contour plot of the effec- 
tive potential in Fig. [2] describe trajectories for scalar 
fields if and (p with different initial values (the col- 



In the case of n scalar fields ipi{i = 1, . . . , n) with 
a flat (fT— model) target space, the slow-roll parameters 
for the spatially flat Friedmann Universe read (see e.g. 

Bill): e = ^EILiCft^^)' ~ k\du\Vu^; m = 

-Cpi/{H^i) = 2dlH/H ^ -e + Y:]^-,af.UdjU/{Ud,U), 

where di := d/dipi and \dU\'^ = E?=i {diUf- In 
some papers (see e.g. [3d |) it was introduced a "cu- 
mulative" parameter rj = — ipiipi/{H\ifi\^) —e + 
j:7j^^{dfjUmU)(djU)/(U\dU\^) , where \<p\^ = ELi 
We can easily find that for the potential l|3.10|l parameter 77 
coincides exactly with parameters 771 and 1^2 ■ 



ored dots). The time evolution of these scalar fields^ 
is drawn in Fig. [31 Here, the time t is measured in the 
Planck times and classical evolution starts a,t t — 1. 
For given initial conditions, scalar fields approach the 
global minimum of the effective potential along spiral 
trajectories. 

We plot in Figure [4] the evolution of the logarithms 
of the scale factor a{t) (left panel) and the evolution of 
the Hubble parameter H{t) (right panel) and in Fig. 
[5] the evolution of the parameter of acceleration q{t). 

Because for initial condition we use the value a{t = 
1) — 1 (in the Planck units), then loga(i) gives the 
number of e-folds: loga{t) = N{t). The Figure |3] 
shows that for considered trajectories we can reach 
the maximum of e-folds of the order of 10. Clearly, 10 
e-folds is not sufficient to solve the horizon and flat- 
ness problems but it can be useful to explain a part of 
the modern CMB data. For example, the Universe in- 
flates by AA^ « 4 during the period that wavelengths 
corresponding to the CMB multipoles 2 < I < 100 
cross the Hubble radius [33|. However, to have the in- 
flation which is long enough for all modes which con- 
tribute to the CMB to leave the horizon, it is usually 
supposed that AA^ > 15 [38l |. 

The Figure|3]for the evolution of the Hubble param- 
eter (right panel) demonstrates that the red, yellow, 
dark blue and pink lines have a plateau H « const. It 
means that the scale factor a{t) has a stage of the De 
Sitter expansion on these plateaus. Clearly, it hap- 
pens because these lines reach the vicinity of the ef- 
fective potential saddle point and spend there some 
time. 



^ We remind that ip describes fluctuations of the internal space 
scale factor and (p reflects the additional degree of freedom of 
the original nonlinear theory. 
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FIG. 3: Dynamical behavior of scalar fields (left panel) and (j> (right panel) with corresponding initial values denoted 
by the colored dots in Fig. [2l 
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FIG. 4: The number of e-folds (left panel) and the Hubble parameter (right panel) for the corresponding trajectories. 



Th e Fig. [S]for the acceleration parameter defined in 
(|A.6p c onfir ms also the above conclusions. According 
to Eq. (jA.Sp . q = 1 for the De Sitter- like behavior. In- 
deed, all these 4 lines have stages qK, \ for the same 
time intervals when H has a plateau. Additionally, 
the magnification of this picture at early times (the 
right panel of the Figure [5]) shows that pink, green 
and blue lines have also a period of time when q is 
approximately consta nt le ss than one: q ~ 0.75. In 
accordance with Eq. (lA.Sp . it means that during this 
time the scale factor a{t) undergoes the power-law in- 
flation a{t) (X with s « 4. This result confirms our 
rude estimates made above for the trajectories which 
go through t he reg ions where the effective potential 
has the form (|3.10p . After stages of the inflation, the 
acceleration parameter starts to oscillate. Averaging q 
over a few periods of oscillations, we obtain q = —0.5. 
Therefore, the scale factor behaves as for the matter 
dominated Universe: a{t) cx t^^^. Clearly, it corre- 



sponds to the times when the trajectories reach the 
vicinity of the effective potential global minimum and 
start to oscillate there. It is worth of noting, that 
there is no need to plot dynamical behavior for the 
equation of state paramete r Lu{t ) because it is linearly 
connected with q (see Eq. (|A.7p ) and its behavior can 
be easily understood from the pictures for q{t). 



As we have seen above for considered quadratic 
model, the maximal number of e-folds is near 10. Can 
we increase this number? To answer this question, we 
shall consider a new model with a higher degree of 
nonlinearity, i.e. the nonlinear quartic model. 
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FIG. 5: The parameter of acceleration (left panel) and its magnification for early times (right panel). There are two 
different form of acceleration with g ~ 1 (De Sitter-like inflation) and q ~ 0.75 (power-law inflation with s « 4) 
accordingly. The averaging of q over a few periods of oscillations results m q = —0.5 which corresponds to the matter 
dominated decelerating Universe. 



IV. NONLINEAR QUARTIC MODEL 

In this section we consider the nonlinear quartic 
model 



f{R) = R + -fR^~ 2Al 



(4.1) 



For this model the scalar field potential (13. 4p reads 
[H: 



U{cb) = -e 



l/3(gA0 „ ^)4/3 _^ 2A, 



Here, the scalar curvature R and scalar field ( 



(4.2) 
' are con- 



nected as follows: e = /' 



1 -f 47i?3 ^ R 



[(e^* - l)/47] 

We are looking for a solution which has a nonneg- 
ative minimum of the effective potential Uef f{f, (f>) 
(|3.ip where potential C/(0) is given by Eq. (|4.2p . If 0o 
corresponds to this minimum, then, as we mentioned 
above (see also [22]), U{(f>()),X and should be pos- 
itive. To get zero minimum of the effective potential, 
these positive values should satisfy the relation of the 
form of C/(0o) = /i = A/2. Additionally, it is 

important to note that positiveness of U{(j)o) results 
in positive expression for R{(l)o) > TTl. 

Eq. (|4.2p shows that potential U{cj)) has the fol- 
lowing asymptotes for positive 7 and A^^^ : (p ^ 



—00 
and 



U{cf>) « ie- 



[1(47) 



1/3 _^ 

(47)-i/3, 



2Ai3] +00 

{-B+4A/3)4, 



00 =^ U{(t>) w |(47j ^'-e^ 
-|-0. For the latter asymptote we took into account 
that -B + 4:A/3 = (£> - 8)/3y/{D - 2){D - 1) < 
for D < 8. Obviously, the total number of dimen- 
sions D = 8 plays the critical role in quartic nonlin- 
ear theories (see [Tl|, O |33|) and investigations for 
D < 8, D = 8 and D > 8 should be performed sepa- 
rately. To make our paper is not too cumbersome, we 
consider the case D < 8 (i.e. di = 1,2,3), postponing 
other cases for our following investigations. 

It is worth of noting that for considered signs of 
parameters, the effective potential Ueff{f, 4>) (|3.ip ac- 
quires negative values when (f) ~> +00 (and U{4>) — > 0). 
For example, if (/{(po) = fi = A/2 (the case of zero 
minimum of the effective potential), the effective po- 
tential Ueffiip.fp ^ 00) < for < e-'"f < 2 
and the lowest negative asymptotic value Ue 
— (16/27)A takes place along the line e~ 
Therefore, zero minimum of Ueff is local^^. 

As we mentioned above, extremum positions (pi of 
the potential [/(0) coincide with extremum positions 



f 'min 

-hip 



4/3. 



= 0. The 



condition of extremum for the potential U{<j)) reads: 



dU 

^ = 



di) 



7(4 -di) 



R + 2Ad 



(4 + rfi) 
7(4 -di) 



(4.3) 



Negative values of A^i and 7 may lead either to negative min- 
ima, resulting in asymptotically AdS Universe, or to infinitely 
large negative values of Uf^ff In the present paper we 

want to avoid both of these possibilities. Therefore, we shall 
consider the case of Ai3,7 > 0. See also footnote 1121 



It is not difficult to show that the thin shell approximation 
is valid for considered model and a tunnelling probability 
from the zero local minimum to this negative U^f f region is 
negligible. 
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For positive 7 and this equation has two real roots: 
^ ' 2 \^ Y {4-di)kVM J 

(4.4) 




where we introduced a dimensionless parameter 

k := 7AI5 , (4.6) 

which is positive for positive 7 and A u , and quantities 
M, uj read 



M 



1 



,1/3 



^ glO/3 (4 + rfl) 

3-2i/3/s(4_di) 



(4.7) 



UJ = fc [-27(4-di)(2 + di)^ 



272(4 - di)2(2 + rfi)4 - 4 • 243fc(16 - d2)3 

(4.8) 

It can be easily seen that for fc > we get a; < and 
M > 0. To have real uj, parameter k should satisfy 
the following condition 



k < 



272(4 -di)2(2 + di)^ 
4-243(16 -^2)3 



(4.9) 



It is not difficult to verify that roots ^0(1,2) are real 
and positive if < fc < fco and they degenerate for 
k ^ ko : -Ro(i.2) (AD/2)^/M. In this limit the 
minimum and maximum of [/(0) merge into an inflec- 
tion point. Now, we should define which of these roots 
corresponds to minimum of U{(j)) and which to local 
maximum. The minimum condition 



d2 [/(</)) 



> =^ 7 [{di + 2) - 47i?3(4 - di)] > 

(4.10) 

results in the following inequality ^2- 

7>0 : {di+2)-A-fRl{A-di)>0. (4.11) 

Thus, the root Rq which corresponds to the minimum 
of U {4>) should satisfy the following condition: 



< i?o < 



/ di 



V47(4-di) 



1/3 



(4.12) 



As we have already mentioned above, the condition U{4>q) > 
leads to the inequality R((po) > 11]. Taking into account 
the condition di < 4, we clearly see that inequality (di + 2) + 
4171^0(4 — di) < for 7 < cannot be realized. This is an 
additional argument in favor of positive sign of 7. 



Numerical analysis shows that i?o(i) satisfies these 
conditions and corresponds to the minimum. For 

i?o(2) we obtain that i?o(2) > ( 47(4-^1) ) corre- 
sponds to the local maximum of U {(p) . In what follows 
we shall use the notations: 



-1- 



47i?i 



47i?i 



and U{(t)rarn) 

note that 



•0(1) 

.3 

•0(2) 
Urai 



(4.13) 
(4.14) 
We should 



1) — Umiji , U {4>rnax) 

^min^ft'max ^nd the ratio U max I U mill de- 
pend on the combination k (|4.6|) rather than on 7 and 
Ad taken separately. 

Obviously, because potential U{4>) has two extrema 
at 0mm and 4>max i the effective potential f/e/ / {ip, (j)) 
may have points of extrema only on the lines (f) = (j^min 
and ip = (t)max where dUef }/ d4>\^^.^^ij,^^^ = 0. To find 
these extrema of Uef / , it is necessary to consider the 
second extremum condition dUef / / dLp = on each 
line separately: 



eff 



dtp 




- 3fhl + 2Axi = , 

- 3/2x1 + 2Ax2 = , 



(4.15) 

where xi = exp (^—^y2di/{di + 2) 1^1^ > and X2 = 

cxp (^—y/2di/{di + 2) > 0; pi and p2 denote 

positions of extrema on the lines (j) = 4'min and 
= 4>max, respectively. These equations have solu- 



tions 



xi(±; 



= a ± \/ o? — (i , 



X2(±) 



a ± 



a > 



u„ 



Urr 



a2 > ai 



(4.16) 



(4.17) 



where we have introduced the notations: a = X/ (3/2) 
and P = t/min/(3/f). These equations show that 
there are 5 different possibilities which are listed in 
the Table HI 

To clarify which of solutions (|4.16p and (|4.17p cor- 
respond to minima of the effective potential (with re- 
spect to p) we should consider the minimum condition 



eff 



d'^p 



> 



f/e.t. + x'9/?-4Ax>0, 



(4.18) 

where Uextr is either Umin or Umax and x denotes 
either xi or X2- Taking into account relations (j4.15p . 
we obtain 



x'S/i^ - xA > 



X > 



A 
3?! 



(4.19) 
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TABLE I: The number of extrema of the effective potential Ueff depending on the relation between parameters. 



< a < ai 


a = ai 


Qi < a < a2 


a = a2 


a > a2 


no extrema 


one extremum 

(point of 
inflection on 
the line 


two extrema 

(one minimum 
and one saddle 
on the line 

<!> = 0min ) 


three extrema 

(minimum and 
saddle on the line 

inflection on 
the line 

= (f>max) 


four extrema 

(minimum and 
saddle on 
the line 

^ — ^mirt 

maximum and 
saddle on 
the line 



Thus, roots Xi.2(+) define the positions of local min- 
ima of the effective potential with respect to the vari- 
able if and Xi,2(-) correspond to local maxima (in the 
direction of if). 

Now, we fix the minimum Xi(+) the point ip = 0. 
It means that in this local minimum the internal space 
scale factor is stabilized at the present day value. In 
this case 



Xi{+)\ 



^0 = 1 = a + y/a^-f^ 



a < 1 



_ 1+/3 
2 

(4.20) 
/3e [0,1). 



Obviously, we can do it only if^^ 
For xi(-) we get: xi(-) = 

Additionally, the local minimum of the effective po- 
tential at the point (ip — 0, (j) = (j)min) should play the 
role of the nonnegative four-dimensional effective cos- 
mological constant. Thus, we arrive at the following 
condition: 



A, 



UeffiV = 0, = <pmin) 
-A + U■^^■a + /i > 

o 



(4.21) 



From the latter inequality and equation (|4.20p we get 
/? € [|,l). It can be easily seen that [3 — l/i (and, 
correspondingly, a = 2/3) results in Ae// = and 
we obtain the mentioned above relations: Umin = 
f1 = A/2. In general, it is possible to demand that 
Kef f coincides with the present day dark energy value 
10~^'^cm~^. However, it leads to very flat local min- 
imum which means the decompactification of the in- 
ternal space [23]. In what follows, we shall mainly 
consider the case of zero Ae// although all obtained 
results are trivially generalized to Ag// = 10~^^cm~^. 

Summarizing our results, in the most interesting 
case of a > a2 the effective potential has four ex- 
trema: local minimum at I (/jI = Q,(j)min], local 



Particular value a = 1 corresponds to the case q = ai = 
where the only extremum is the inflection point with Xi(-) ■ 
Xi{+) = o = 1. Here, A = Umm = 3/f and Ueff (if = 0,4>- 

4'rriin) = — A 4- Umin + /l > 0- 



maximum at (99 1 



X2(-) J vmax) and two saddle-points at 
(see Fig. [T]). 



and 



X2( + ) 



We pay particular attention to the case of zero local 
minimum J7e//(<^|^^^^^ = 0,0™™) = where [3 = 

1/3 =^ a (1 + /3)/2 = 2/3. To satisfy the four- 
extremum condition a > q;2, we should demand 



4 
<3 



(4.22) 



The fraction Umax/Umm is the function of k and 
depends parametrically only on the internal space 
dimension c?i. Inequality (j4.22p provides the lower 
bound on k and numerical analysis (see Fig. [S]) gives 

Umax/Umin 
3 




0.005 0.01 0.015 0.02 



0.5- 



FIG. 6: The form of Umax/Umin as a function of fc G (0, fco] 
for di = 1, 2, 3 from left to right, respectively. Dashed line 
corresponds to Umax/Umin = 4/3. 

k{di = 1) w 0.000625 ; k{di = 2) w 0.00207 ; k{di = 
3) « 0.0035 . Therefore, effective potentials with zero 

local minimum will have four extrema if k G (fc, fco) 
(where fcp is defined by Eq. (|4.9p ). The limit k ^ k 
results in merging X2(-) ^ X2(+) and the limit k ^ ko 
resuhs in merging xi(-) ^ X2(-) and xi(+) ^ X2(+)- 
Such merging results in transformation of correspond- 
ing extrema into inflection points. For example, from 
Fig. [5] follows that Umax/Umm ^ 1 for fc ^ fco. 
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The typical contour plot of the effective potential 
with four extrema in the case of zero local minimum 
is drawn in Fig. [71 Here, for di = 3 we take k — 
0.004 e (fc, ko) which gives a2 « 0.655. Thus, a = 
2/3 « 0.666 > a2. 

Let us investigate now a possibility of inflation for 
considered potential. First of all, taking into account 
the comments in previous section (see a paragraph be- 
fore Eq. (|3.9p ). it is clear that topological inflation in 
the saddle point Xi(-) as well as the slow rolling from 
there in the direction of the local minimum Xi(+) are 
absent. It is not difficult to verified that the gener- 
alized power-low inflation discussed in the case of the 
nonlinear quadratic model is also absent here. Indeed, 
from Eqs. p.ip and (|4.2p follows that nonlinear po- 
tential U{(j>) can play the leading role in the region 
(j) — oo (because U{(j)) — > for (/) ^ -l-oo). In this 
region Ueff oc exp (— y^(^) exp (— .^^0) where q = 
2di/idi+2) and q2 = = [di+Af /[{di + 2){di+?,)]. 
For these values of q and q2 the slow-roll conditions 
are not satisfied: e « 771 « 772 ~ g/2 -1-92/2 > 1. How- 
ever, there are two promising regions where the stage 
of inflation with subsequent stable compactification of 
the internal space may take place. We mean the local 
maximum X2(-) f^nd the saddle X2(-i-) (see Fig. [7]). Let 
us estimate the slow roll parameters for these regions. 

We consider first the local maximum X2(-)- It is ob- 
vious that the parameter e is equal to zero here. Addi- 
tionally, from the form of the effective potential (|3.ip 
it is clear that the mixed second derivatives are also 
absent in extremum points. Thus, the slow roll param- 
eters 771 and 772, defined in the footnote ([5]), coincide 
exactly with rj^ and 77^. In Fig. [5] we present typical 

form of these parameters as functions of fc e (A:, fcg) in 
the case (3 = 1/3 and di = 1,2,3. These plots show 
that, for considered parameters, the slow roll inflation 
in this region is possible for di — 1,3. 

The vicinity of the saddle point X2(+) is another 
promising region. Obviously, if we start from this 
point, a test particle will roll mainly along direction 
of (j). That is why it makes sense to draw only \rirj,\. In 
Fig. [9] we plot typical form of 1 7^0! in the case (3 — 1/Z 
and di = 1, 2, 3. Left panel represents general behav- 
ior for the whole range of /c € {k,ko) and right panel 
shows detailed behavior in the most interesting region 
of small k. It shows that di = 3 is the most promising 
case in this region. 

Now, we investigate numerically the dynamical be- 
havior of scalar fields and the external space scale 
factor for trajectories which start from the regions 
Xi(-) , X2(-) and X2(+) • AH numerical calculations per- 
form for /? = 1/3, di = 3 and A: = 0.004. The colored 
lines on the contour plot of the effective potential in 
Fig. [7| describe trajectories for scalar fields (p and (p 
with different initial values (the colored dots) in the 
vicinity of these extrema points. The time evolution 
of these scalar fields is drawn in Fig. [TOl For given 
initial conditions, scalar fields approach the local min- 
imum Xi(+) of the effective potential along the spiral 
trajectories. 

We plot in Figure [TT] the evolution of the logarithm 



of the scale factor a{t) (left panel) which gives directly 
the number of e-folds and the evolution of the Hub- 
ble parameter H{t) (right panel) and in Fig. [T^ the 
evolution of the parameter of acceleration q{t). 

The FigurefTTjshows that for considered trajectories 
we can reach the maximum of e-folds of the order of 22 
which is long enough for all modes which contribute 
to the CMB to leave the horizon. 

The Figure [11] for the evolution of the Hubble pa- 
rameter (right panel) demonstrates that all lines have 
plateaus H « const. However, the red, yellow and 
blue lines which pass in the vicinity of the saddle X2{+) 
have bigger value of the Hubble parameter with re- 
spect to the dark blue line which starts from the Xi(-) 
region. Therefore, the scale factor a{t) has stages of 
the De Sitter-like expansion corresponding to these 
plateaus which last approximately from 100 (dark blue 
line) up to 800 (red line) Planck times. 

The Fig. [T2| for the acceleration parameter con- 
firms also the above conclusions. All 4 lines have 
stages g « 1 for the same time intervals when H has 
plateaus. After stages of infiation, the acceleration 
parameter starts to oscillate. Averaging q over a few 
periods of oscillations, we obtain q = —0.5. Therefore, 
the scale factor behaves as for the matter dominated 
Universe: a{t) cx t^/"^. Clearly, it corresponds to the 
times when the trajectories reach the vicinity of the 
effective potential local minimum Xi(-i-) ^^^d start to 
oscillate there. 

Let us investi gate now a possibility of the topo- 
logical inflation [161 [40j | if scalar fields ip, 4> stay in 
the vicinity of the saddle point X2(+)- As we men- 
tioned in Section 2, topological inflation in the case 
of the double-well potential takes place if the dis- 
tance between a minimum and local maximum big- 
ger than A0cr = 1.65. In this case domain wall 
is thick enough in comparison with the Hubble ra- 
dius. The critical ratio of the characteristic thickness 
of the wall to the horizon scale in local maximum is 
r^,H w |J7/3a[/00|i/2 w 0.48 [13 and for topological 
inflation it is necessary to exceed this critical value. 
Therefore, we should cheque the saddle X2(+) from 
the point of these criteria. 

In Fig. [131 (left panel) we draw the difference 
= (jiraax " (prmn for the profile Lp = v\x2{+) as a 
functions of fc G {k,ko) in the case (3 — 1/3 for di- 
mensions di = 1,2,3. This picture shows that this 
difference can exceed the critical value if the num- 
ber of the internal dimensions is di = 2 and di = 3. 
Right panel of Fig. [T3| confirms this conclusion. Here 
we consider the case /3 = 1/3, k = 0.004 and di = 3. 
For chosen values of the parameters, A0 = 2.63 which 
is considerably bigger than the critical value 1.65 and 
the ratio of the thickness of the wall to the horizon 
scale is 1.30 which again bigger than the critical value 
0.48. Therefore, topological inflation can happen for 
considered model. Moreover, due to quantum fluctu- 
ations of scalar fields, inflating domain wall will have 
fractal structure: it will contain many other inflating 
domain walls and each of these domain walls again 
will contain new inflating domain walls and so on [la | . 
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FIG. 7: Contour plot of the effective potential Uef f{f, <j)) P-ip with potential (7(0) of the form (|4.2p for parameters 
(3 = 1/3, di = 3 and k — 0.004. This plot shows the local zero minimum, local maximum and two saddles. The colored 
lines describe trajectories for scalar fields starting at different initial conditions. 



^^1 1^^ 




FIG. 8: Graphs of |77^| (left panel) and \ri^\ (right panel) as functions of fc G (fc,fco) for local maximum X2(-) and 
parameters /3 = 1/3 and di = 1, 2, 3. 



Thus, from this point, such topological inflation is the 
eternal one. 

To conclude this section, we want to draw the at- 
tention to one interesting feature of the given model. 
From above consideration follows that in the case of 
zero minimum of the effective potential the positions 
of extrema are fully determined by the parameters k 
and di, and for fixed k and di do not depend on the 
choice of Ad. The same takes place for the slow roll 
parameters. On the other hand, if we keep k and di, 
the hight of the effective potential is defined by 
(see Appendix B). Therefore, we can change the hight 
of extrema with the help of A^i but preserve the con- 
ditions of infiation for given k and c?i . 

However, the dynamical characteristics of the model 
(drawn in figures [TU1- [T^ depend on variations of A^i 
by the self-similar manner. It means that the change 



of hight of the effective potential via transformation 
Ao — !■ cAd (c is a constant) with fixed k and di results 
in rescaling of figures (TU] - HH in l/\/c times along the 
time axis. 



V. SUMMARY AND DISCUSSION 

In our paper we investigated a possibility of infla- 
tion in multidimensional cosmological models. The 
main attention was paid to nonlinear (in scalar curva- 
ture) models with quadratic and quartic la- 
grangians. These models contain two scalar fields. 
One of them corresponds to the scale factor of the in- 
ternal space and another one is related with the non- 
linearity of the original models. The effective four- 
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FIG. 9: Grafs of \ri^\ as functions of k for saddle point X2{+) and parameters /3 = 1/3 and di — 1,2,3. Left panel 
demonstrates the whole region of variable k £ (k, ko) and right panel shows detailed behavior for small k. 




FIG. 10; Dynamical behavior of scalar fields (f (left panel) and (right panel) with corresponding initial values denoted 
by the colored dots in Fig. [T] 
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FIG. 11: The number of e-folds (left panel) and the Hubble parameter (right panel) for the corresponding trajectories. 
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FIG. 12: The parameter of acceleration (left panel) and its magnification for early times (right panel). 



0max-0 



mm 



U 




eff 



0.005 0.01 0.015 0.02 




FIG. 13: Left panel demonstrates the difference (pmax — (jtmin (for the profile ip — 'p\x2(+)) ^ ^ functions of fc G (fc,fco) 
for parameters /3 — 1/3, and di — 1,2,3 (from left to right respectively). Dashed line corresponds to (pmax — <f>min = 
1.65. Right panel shows the comparison of the potential L''e//(v5|x2(+) i with a double-well potential for parameters 
/3 = 1/3, fc = 0.004 and di = 3. 



dimensional potentials in these models are fully de- 
termined by the geometry and matter content of the 
models. The geometry is defined by the direct prod- 
uct of the Ricci-flat external and internal spaces. As 
a matter source, we include a monopole form field, D- 
dimensional bare cosmological constant and tensions 
of branes located in fixed points. The exact form of 
the effective potentials depends on the relation be- 
tween parameters of the models and can take rather 
complicated view with a number of extrema points. 

First of all, we found a range of parameters which 
insures the existence of zero minima of the effective 
potentials. These minima provide sufficient condition 
to stabilize the internal space and, consequently, to 
avoid the problem of the fundamental constant vari- 
ation. Zero minima correspond to the zero effective 
four-dimensional cosmological constant. In general. 



we can also consider positive effective cosmological 
constant which corresponds to the observable now 
dark energy. However, it usually requires extreme fine 
timing of parameters of models. 

Then, for corresponding effective potentials, we in- 
vestigated the possibility of the external space infla- 
tion. We have shown that for some initial conditions 
in the quadratic and quartic models we can achieve 
up to 10 and 22 e-folds, respectively. An additionally 
bonus of the considered model is that model can 
provide conditions for the eternal topological inflation. 

Obviously, 10 and 22 e-folds are not sufficient to 
solve the homogeneity and isotropy problem but big 
enough to explain the recent CMB data. To have 
the infiation which is long enough for modes which 
contribute to the CMB, it is usually supposed that 
AA^ > 15 [3^. Moreover, 22 e-folds is rather big 
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number to encourage the following investigations of 
the nonlinear multidimensional models to find theories 
where this number will approach 50-60. We have seen 
that the increase of the nonlinearity (from quadratic 
to quartic one) results in the increase of AA'^ in more 
that two times. So, there is a hope that more compli- 
cated nonlinear models can provide necessary 50-60 
e-folds. Besides, this number is reduced in models 
where long matter dominated (MD) stage followed in- 
flation can subsequently decay into radiation [23, 123 • 
Precisely this scenario takes place for our models. 
We have shown for quadratic and quartic nonlinear 
models, that MD stage with the external scale factor 
a ~ t^^"^ takes place after the stage of inflation. It 
happens when scalar fields start to oscillate near the 
position of zero minimum of the effective potential. 
However, scalar fields are not stable. For example, 
scalar field ip decays into two photons (p —^ 2^ with 



the decay rate F 



m: 



V 



{Mpi/m^ftpi. 



30( 1 . Thus the life time 
The reheating tempera- 



ture is given by the expression Tpn ^ (w^ A^^p;)^^^- 
Therefore, to get Tpn ^ IMeV necessary for the nu- 
cleosynthesis, we should take m^p > lOTeV. In paper 
[43j |. it is shown that for such scenario with interme- 
diate MD stage, the necessary number of e-folds is 
reduced according to the formula: 



1 
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AN = -gln(^y.97 
1 /45 _ 



-3/2 


ml 




TMpi 


3/2 


Mpi\ 




J 



(5.1) 
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APPENDIX A: FRIEDMANN EQUATIONS 
FOR MULTI-COMPONENT SCALAR FIELD 
MODEL 



We consider n scalar fields minimally coupled to 
gravity in four dimensions. The effective action of 
this model reads 



1 



167rG 



;;(o)l 



- G^,fg^''>'^''d^ip^d,ip^ ~ 2U{v\^\ . . .)) (A.l) 

where the kinetic term is usually taken in the canon- 
ical form: dj = diag(l, 1, . . .) (flat a model). Such 
multi-component scalar fields originate naturally in 
multidimensional cosmological models (with linear or 
nonhnear gravitational actions) [3, [3,[T3]. We use the 
usual conventions c = ?i=l, i.e. Lpi = tpi = 1/Mpi 
and SttG = Sn/Mpi. Here, scalar fields are dimension- 
less and potential U has dimension [U] — lehgth^^. 

Because we want to investigate dynamical behav- 
ior of our Universe in the presence of scalar fields, 
we suppose that scalar fields are homogeneous: (p* = 
(p^{t) and four-dimensional metric is spatially- flat 
Friedmann-Robertson- Walker one: 5^"-' ~ ~dt ® dt + 
a?{t)dx (g) dx. 

For energy density and pressure we easily get: 



where counts the effective number of relativistic 
degrees of freedom and we took into account that de- 
caying particles are scalars. This expression weakly 
depends on g^. For example, if ~ lOTeV we ob- 
tain -6.27 < A7V < -5.11 for 1 < < 10^. Thus, 
AA^ « —6. Therefore, we believe that the number 
of e-folds is not a big problem for multidimensional 
nonlinear models. The main problem consists in the 
spectral index. For example, in the case of i?"* model 
we get Ug 1 + 2ry|^2(+) ~ 0-61 which is less than 
observable now « 1. A possible solution of this 
problem may consist in more general form of the non- 
linearity /(-R). It was observed in [l^ that simultane- 
ous consideration quadratic and quartic nonlinearities 
can flatten the effective potential and increase . We 
postpone this problem for our following investigations. 
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SttG \2 

1 ^l 



(A.2) 
(A.3) 



SttG \2 

( \G,,^'iP^ ^A^G{p + P) 
\u = AirGip - P) . 
The Friedmann equations for considered model are 



3 ( ^ j =3H^ = SnGp = ^G^j^V^ + U , (A.l) 



and 
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H = -AttG{p + P) = --G,,^V-'' 



(A.5) 



From these 2 equations, we obtain the following ex- 
pression for the acceleration parameter: 

a A-kG , , SttG, 

1 ^ 7777 = l-^(P + ^) = -^(/' + 3P) 



ma 
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MP 



-4 X iGy^V' + 2C/ 
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(A.6) 
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It can be easily seen that the equation of state (EoS) 
parameter lo = P/ p and parameter q are linearly con- 
nected: 



:(l + 3u;). 



(A.7) 



From the definition of the acceleration parameter, it 
follows that q is constant in the case of the power-law 
and De Sitter-like behavior: 



_ r (,s-i)/s 



1 



a cx e 



a cx , 
Ht 



(A., 



For example, q — —0.5 during the matter dominated 
(MD) stage where s = 2/3. 

Because the minisuperspace metric Gij is flat, the 
scalar field equations are: 



dU 



0. 



(A.9) 



For the action (jA.ip . the corresponding Hamilto- 
nian is 



where 



Pr 



^^^^^ 



(A.IO) 



(A.ll) 



are the canonical momenta and equations of motion 
have also the canonical form 

It can be easily seen th at the latter equation (for Pi) 
is equivalent to the eq. (jA.9[) . 

Thus, the Friedmann equations together with the 
scalar field equations can be replaced by the system 
of the first order ODEs: 
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H = 



a 

i (-4 X iCy-^V-'' +2U]-H^ (A.16) 



with Eq. (|A.4[) considered in the form of the initial 
conditions: 



H{t = 0) 
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(A.17) 



We can make these equations dimensionless: 
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dip' 
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(A.18) 



(A.19) 



That is to say the time t is measured in the Planck 
times tpi , the scale factor a is measured in the Planck 
lengths Lpi and the potential U is measured in the 
Afp; units. 

We use this system of dimensionless f irst o rder 
ODEs together with the initial condition (|A.17p for 
numerical calculation of the dynamics of considered 
models with the help of a Mathematica package [361] . 



APPENDIX B: SELF-SIMILARITY 
CONDITION 



Due to the zero minimum conditions C/(0min) = 
/2 ^ A/2, the effective potential ((XT|) can be written 
in the form: 



U{(j)min)e 



rjdT 



' V + 



2e 



nidi 
V + 



.(B.l) 



Exact expressions for C/(0) p.7p and (|4.2p indicate 
that the ratio 



F{(j),k,di) 



(B.2) 



depends only on (f>, k and di. Dimensionless parameter 
k = S^Ajj for the quadr atic m odel and k = for the 
quartic model. In Eq. (|B.2[) we take into account that 
(prnin is a function of k and di: (j^min = 0mm(fc,(ii)- 
Then, C/(0m„) defined in Eqs. (1^ and (lO)) reads: 



ADF{(f)rmn{k,di),k,di) . 



(B.3) 



t=0 



Therefore, parameters k and di determine fully the 
shape of the effective potential, and parameter A^ 
serves for conformal transformation of this shape. 
This conclusion is confirmed also in sections 3 and 
4 where we show that positions of all extrema de- 
pend only on k and di. Thus, figures [2l and [7] for 
contour plots are defined by k and di and will not 
change with A^i. From the definition of the slow roll 
parameters it is clear that they also do not depend 
on the hight of potentials and in our model depend 
only on k and di (see figures [8] and |9]) . Similar depen- 
dence takes place for difference A(f> — (pmax — 4>min 
drawn in Fig. 1131 Thus the conclusions concerning 
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the slow roll and topological inflations arc fully deter- 
mined by the choice of k and c?i and do not depend on 
the hight of the effective potential, in other words, on 
A/). So, for fixed k and di parameter A^i can be arbi- 
trary. For example, we can take in such a way that 
the hight of the saddle point X2{+) will correspond to 
the restriction on the slow roll inflation potential (see 
e.g. [SI) C/e// < 2.2 X 10~"A/|;, or in our notations 
Ueff < 5.5 X 10-^°M^i. 

Above, we indicate figures which (for given k and 
di) do not depend on the hight of the effective po- 
tential (on Ad). What will happen with dynamical 
characteristics drawn in figures [TOl E] and [T^] (and 
analogous ones for the quadratic model) if we, keep- 
ing fixed k and di, will change A^i? In other words, 
we keep the positions of the extrema points (in (tp, (p)- 
plane) but change the hight of extrema. We can easily 
answer this question using the self-similarity condi- 



tion of t he Fr iedm ann e quations. Let the potential U 
in Eqs. (jA.2[) and (jA.Sp be transformed conformally: 
U ^ cU where c is a constant. Next, we can intr oduc e 
a new time variable r := ^/ct. Then, from Eqs. (jA.2[) - 
(jA.5|) follows that the Friedmann equations have the 
same form as for the model with potential U where 
time t is replaced by time r. This condition we call 
the self-similarity. Thus, if in our model we change 
the parameter A^ : A^i — > cA^, it results (for fixed 
k and di) in rescaling of all dynamical graphics (e.g. 
Figures \W\ - [T2|) along the time axis in 1/ y/c times 
(the decrease oi Ad leads to the stretch of these fig- 
ures along the time axis and vice versa the increase of 
Ad results in the shrink of these graphics). Numeri- 
cal calculations confirm this conclusion. The property 
of the conformal transformation of the shape of Uef / 
with change of A^i for fixed k and di can be also called 
as the self-similarity. 
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